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ABSTRACT 

Current-voltage characteristics have been obtained from Langmuir 
probes at various positions inside a Kaufman-type ion motor. Data are 
presented for the engine in static (no beam, but plasma) and operat- 
ing (160-ma beam) conditions. The analysis of the probe data 
indicates that the electrons in the plasma possess a non-Maxwellian 
distribution of energies. To obtain plasma properties from the probe 
curves, the non-Maxwellian energy distribution is assumed to consist 
of a primary monoenergetic group of high-energy electrons plus the 
usual Maxwell distribution. Working from this assumption the probe 
curves are analyzed for primary and Maxwell energies and densities, 
and the data for radial positions inside the engine are presented. It is 
shown from the density curves that the primary electrons are the prin- 
cipal source of ion production. The beam current as calculated from 
probe data near the exit grids agrees to the correct order of magnitude 
with the measured value. 


I. INTRODUCTION 


The possibilities of the Kaufman-type electron bom- 
bardment ion source ( Ref. 1 ) for use as a thrust device 
for deep-space applications is well documented ( Ref. 2 ) . 
Even though the engine in its present state achieves 
reasonable efficiencies, its operation is not thoroughly 
understood. It has been studied externally; the variation 
of beam current, propellant utilization efficiency, etc., 
have been measured as the various external operating 
parameters such as arc voltage and magnetic field have 
been changed. For a deeper understanding of the source 
operation, it is necessary to probe the interior of the 
engine, where the. plasma is being formed. Only in this 
way can a theoretical model be correlated with the 
actual situation. Then, if a proper model can be found, a 
calculation of the efficiencies and beam properties can 
be made from the operating parameters, and it will be 
possible to attempt to improve the engine operation for 
use as a thrust device. 


This paper deals with an experimental program set up 
to investigate the properties of the plasma inside an 
operating Kaufman-type engine. From the results of this 
investigation it is hoped that a theoretical model can be 
developed which will sufficiently explain plasma forma- 
tion so that future engines of this type can be built with 
predictable characteristics. 


Owing to their small size and ability to distinguish 
plasma properties at different points, Langmuir probes 
(Ref. 3) were used to probe the inside of the engine. 
These probes are simple to operate and need a minimum 
amount of circuitry, advantages not to be overlooked 
when the engine is operated under flight-simulated con- 
ditions during which the whole device is floating at a 
potential of about 5 kv. 
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II. EXPERIMENTAL SETUP 


For this investigation an ion engine of the Kaufman 
type was purchased from Goodrich-High Voltage Astro- 
nautics, Inc, 1 A schematic diagram of the engine and 
associated power supplies is given in Fig. 1. 

The propellant, mercury, enters the arc or discharge 
chamber as a gas through a distributor in the rear of the 
engine. A cylindrical, nickel-matrix oxide, indirectly 
heated cathode is positioned along the axis of the engine. 
A cylindrical anode surrounds the cathode and is attached 
to, but electrically isolated from, the engine housing. 
The accelerator electrodes consist of two 0.060-in. molyb- 
denum plates through which a matrix of 0.187-in. holes 
has been drilled in a close-pitch hexagonal array. One 
of the electrodes is attached to and is electrically com- 


J Now Ion Physics Corp., a subsidiary of High Voltage Engineering 
Corp., Burlington, Mass. 


mon with the engine housing which forms the arc cham- 
ber. The second electrode is mounted downstream of the 
first on three ceramic insulators. A solenoidal magnet is 
wound around the entire housing to provide field diver- 
gence downstream of the accelerating electrodes. A 
decelerating electrode is not employed. 

During engine operation, the housing, cathode sheath, 
magnet, and propellant feed system are electrically com- 
mon and maintained at high positive potential, while the 
outer accelerating electrode is biased negatively, typi- 
cally — 1000 to —1500 v. The anode remains positive 
with respect to engine housing by the arc or discharge 
potential. The cathode temperature is varied to achieve 
a desired arc current for any selected arc voltage, thus 
ensuring emission-limited cathode operation. 

The engine was installed in a 3- by 7-ft vacuum cham- 
ber capable of being evacuated to 10" 7 mm Hg (Fig. 2). 



SYSTEM 

GROUND 


Fig. 1. Circuit diagram of engine and probe 
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A cryogenically cooled collector was positioned 30 in. 
downstream of the engine-accelerating electrode. The 
collector is a two-section type, providing beam current 
measurements and suppression of secondary electrons. 



Fig. 2. View of experimental setup 


A movable Langmuir probe (Fig. 3) was constructed 
by encapsulating a 0.020-in. -diameter pure tungsten wire 
within a 96% fused silica sheath packed with finely 
powdered boron nitride. The tungsten wire was allowed 
to protrude 0.300 in. from the end of the sheath. 

The probe was inserted into the plasma from the rear 
of the engine through a boron nitride seal and connected 
to a traversing and radial positioning mechanism located 
outside the vacuum tank. This mechanism was designed 
and constructed to provide probe position information 
to a tolerance of ±0.015 in. in both the axial and radial 
directions. 

Probe instrumentation, which consisted of a 0- to 250- 
vdc power supply, variable frequency function generator, 
cathode follower amplifier, x-y plotter with an input 
impedance of 1 megohm, oscilloscope, bias batteries, and 
100-ohm shunt, was electrically connected together as 
illustrated in Fig. 1. The instrumentation was isolated 
from ground potential and coupled to the ion engine and 
probe through two high-voltage vacuum relays, thus 
providing instrumentation decoupling capability without 
disturbing engine operation. All probe voltage infor- 
mation was referenced to anode potential. 


A second Langmuir probe of variable area (Fig. 4) 
was constructed, for reasons discussed later in this paper. 
This probe consisted of a 0.090-in. -diameter pure tung- 
sten rod sheathed in a boron nitride insulator. One end 
of the tungsten rod was attached to a solenoid which, 
when actuated, could vary the length of the rod protrud- 
ing from the insulator. In the retracted position the rod 
protruded 0.050 in. from the insulator; in the extended 
position it protruded 0.350 in. The probe was inserted 
into the plasma from the rear of the engine. No travers- 
ing mechanism was employed on the second probe. 



Fig. 3. Movable probe. Insert shows approximate probe 
positioning points inside engine 



Fig. 4. Variable-area probe 
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III. TYPICAL RESULTS 


The analysis of a Langmuir probe characteristic (LPC) 
is simple if the electrons in the plasma are distributed 
according to the Maxwell-Boltzmann law. In this case a 
plot of the LPC on semilogarithmic graph paper will 
yield a straight line. 

Quite early in this program many LPCs were obtained 
which did not plot as straight lines. Sometimes the 
deviation was small, sometimes large. In most cases the 
upper portion of the LPC was linear, the deviation from 
linearity not occurring until the probe was some 5 to 10 v 
below plasma potential. In such cases, the question arises 
as to whether the deviations are due to an actual non- 
Maxwellian energy distribution law, or to some external 
cause such as an improper circuitry or a probe surface 
contamination, or to one of the many other phenomena 
evoked by experimenters using Langmuir probes. A 
definite answer to the question is hard to give, but from 
a series of independent experiments (some to be de- 
scribed herein ) it appears that the anomalous LPCs are a 
result of a non-Maxwellian distribution law for the 
electrons in the plasma. It will also become apparent that 
the deviation from Maxwellian may be quite important 
in determining the operation of the source. A similar 
non-Maxwellian distribution has been reported by other 
workers, notably Medicus (Ref. 4), who studied a fire- 
ball discharge in neon. 

Figure 5 presents a typical semilogarithmic plot of an 
LPC taken with the movable probe. Note the deviation 
from a straight line at negative voltage values. This type 
of behavior is typical of most of the LPCs taken under 
various operating conditions inside the engine. However, 
if the arc voltage is low enough, as shown in Fig. 5, the 
LPC follows a straight line for three decades of the 
current. 

It is known that for some plasmas an anomalous 
characteristic such as that shown in Fig. 5 is found when 
the plasma creeps up the insulator of the probe with 
increasing negative voltages. To show that the anomalous 
LPCs do not occur from some probe-insulator-plasma 
interaction, the variable-area probe described previously 
was used. Figure 6 illustrates the results. The curve 
labeled ‘extended” is the LPC obtained by the probe in 
the extended position; any probe-insulator interaction 
should be present in this curve. The curve labeled 
“difference” is the difference between the two LPCs 


obtained by the probe in the extended and retracted 
positions; any insulator effect should not appear in_ the 
“difference” curve. A comparison of the “extended” and 
“difference” curves shows no difference, thus ruling out 
any effect of insulator-plasma interaction on the probe 
traces. 



PROSE POTENTIAL, v 


Fig. 5. Typical semilogarithmic LPC showing variation in 
probe trace with arc voltage 

arc voltage— noted magnetic field— 20 gauss (6 amp) 

arc current— 2 amp flowrate— 1.7 gm/hr 

To show that the anomalous LPCs do not result from 
the formation of oxides or amalgams or other contamina- 
tion of the probe surface, the voltage on the probe was 
driven to high positive values. By sweeping the voltage 
on the probe at 0.01 cps, electron currents large enough 
to heat the probe to incandescence could be maintained 
for some tens of seconds. Heating by electron bombard- 
ment should raise the surface temperature of the probe 
high enough to boil off any contaminants. 

It was found that LPCs taken with a probe immedi- 
ately after cleaning in the above manner were identical 
to those taken when the probe was not heated to incan- 
descence. The variable-area probe was not heated to 
incandescence. 

From the above discussion, it is seen that the anoma- 
lous LPCs are probably due to an actual non-Maxwellian 
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energy distribution for the electrons. Since the LPC at 
any voltage V is an integral of the energy distribution 
taken from the energy eV to infinity, an excess of high- 
energy electrons in the distribution would yield an LPC 
of the type observed. This interpretation is made through- 
out the remainder of this paper. Quantitative evaluation 
of the probe curves is made in Section IV. 



-15 -10 -5 0 5 10 


PROBE POTENTIAL, v 

Fig. 6. Typical semilogarithmic LPC 
from variable-area probe 

arc voltage — 20 v magnetic f iel cJ — 1 0 gauss (3 amp) 

arc current— 0.85 amp flowrate— 1.7 gm/hr 

Figure 7 shows five LPCs taken at different radial 
positions in the engine. (Engine operating parameters 
are noted in the caption.) Observe the gradual decrease 
in electron saturation current and general straightening 
of the curves near the anode. The saturation current is 
an approximate measure of the density of electrons 
which, as expected, decreases radially. This type of 
behavior is typical of all the LPCs obtained under many 
different operating conditions and was observed very 
near the grids when accelerating voltage was on or off. 



-30 -20 -10 0 10 20 

PROBE POTENTIAL, v 

Fig. 7. Typical radial semilogarithmic LPCs 


arc voltage— 50 v magnetic field— 14 gauss (4 amp) 

arc current-2 amp flowrate-1 .7 gm/hr 

Under some operating conditions severe Tiash” in the 
saturated electron current portion of the LPCs was ob- 
served when the characteristic was viewed on an oscillo- 
scope. This phenomenon, which was interpreted as some 
sort of plasma oscillation, could be observed by coupling 
the oscilloscope to any portion of the engine. Figures 
8 and 9 illustrate the magnitude of these oscillations. 
The only difference between the two pictures is the 
magnetic field intensity. Figure 8 was taken at 20 gauss; 
Fig. 9 at 25 gauss. These strong oscillations were always 
present above a certain critical field strength and may 
play an important role in engine operation. Since the 
LPCs to be used for analysis were taken with an x-y 
recorder, no meaningful traces could be obtained when 
the strong oscillations were present. Therefore, it was 
not possible to obtain any information on the plasma 
when the magnetic field was above about 20 gauss. 
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Fig. 8. Typical LPC obtained on oscilloscope 

arc voltage— 40 v magnetic field-20 gauss (6 amp) 

arc current— 2 amp flowrate— 1.7 gm/hr 



Fig. 9. Typical LPC obtained on oscilloscope 
showing oscillations 

arc voltage— 40 v magnetic field— 25 gauss 19 amp) 

arc current-2 amp flowrate-1.7 gm/hr 
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IV. ANALYSIS OF LPC 


In order to obtain useful information on the plasma, a 
method of breaking down the LPC into an energy dis- 
tribution for the electrons was devised. Before detailing 
this method, some general criteria are presented showing 
the applicability of Langmuir probes to investigations of 
plasma of the kind found in the Kaufman engine. 

1. For the probe to act like an infinite planar probe, 
the sheath thickness should be small in comparison 
to the size of the probe. For a current density of 
0.004 amp/cm 2 and a voltage of 21.5 v, the sheath 
will be 9.8 by 10 -3 cm, as calculated by the space 
charge equation, which is much smaller than the 
diameter of any of the cylindrical probes used. 

2. Elementary probe theory assumes that there are 
very few collisions in the sheath. This means that 
the sheath must be small in comparison to the mean 
free path of the electrons. In the Kaufmann engine, 
particle densities are of the order of 10 11 particles/ 
cc, which gives a mean free path of several meters 
— much larger than the sheath indeed. 

3. Since the probe is collecting electrons in the mag- 
netic field of the engine, the effect of the magnetic 
field upon collection must either be corrected for 
or must be small enough to be neglected. If the 
cyclotron radius of the electrons is large compared 
to the probe dimensions, the effect of the magnetic 
field will be small. The cyclotron radius for a 5-ev 
electron in a 20-gauss field is 0.3 cm, which is large 
compared to the probe diameter of 0.05 cm. 

From the above data it is clear that the usual infinite 
planar probe theory will be applicable. 

When the plasma electrons possess a Maxwell distri- 
bution of energies, the electron current to an infinite 
planar probe for voltages negative with respect to the 
plasma potential is: 



where 

A = area of probe 
e ~ electronic charge 
N m = number of Maxwell electrons 
k — Boltzmann’s constant 
T e — electron temperature 


M e = mass of the electron 
V = V* - V 0 
V 8 — plasma potential 
V 0 = potential applied to probe 

When the probe is at plasma potential, V = 0 and the 
saturation electron current tsm is 

* - = AeN m ylSfc 

From (1), the current-voltage curve will plot as a straight 
line on semilogarithmic graph paper, and the electron 
temperature can be determined from the slope of this 
line. The plasma potential can be found from the break- 
point in the current-voltage curve. Knowing the plasma 
potential and the electron temperature, the density of 
the Maxwell electrons can be found using (2). 


Now consider the case of a monoenergetic group of 
electrons of number N p and energy E p . Assume the veloc- 
ities to possess a random distribution in space. The elec- 
tron current to an infinite planar probe now is (Ref. 5): 



where V p is the electron-volt energy equivalent of E p . 
Thus, in this case, the current-voltage probe curves will 
be linear. Again, at plasma potential, when V = 0: 




AeN p f2eV~ p 

4 V M e 


(4) 


and since V p can be determined from the current-voltage 
curve, N p can be found from (4). Now, the total probe 
current i, when the two groups of electrons are present, 
is; 


i = ip + i M 

Since V p > kT e /e when the arc voltage is large, for large 
negative voltages V 

I'M ^ ^ Ip 

so that i is approximately linear. This linearity gives an 
indication of how to begin to apply these results to the 
probe curves. 

Figure 10 shows a typical anomalous LPC on a linear 
scale. Note that, for large negative voltages, the curve 
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can be approximated by a straight line corresponding to 
monoenergetic electrons of energy of approximately 22 v. 
When the straight line indicated is subtracted from the 
original data, the remaining curve plots as a straight line 
on a semilogarithmic plot. 2 Figure 11 shows the original 
data on a semilogarithmic plot together with the remain- 
ing curve after the straight line has been subtracted. It is 
seen that the electron energy distribution leading to the 
anomalous LPC may be approximated by a distribution 
of electrons composed of two parts: 

1. A “primary” distribution of monoenergetic elec- 
trons of high energy with randomly oriented 
velocities. 

2. A Maxwell distribution of electrons with a low 
mean energy. 



Fig. 10. Typical LPC from x-y recorder (linear scale) 

ore voltage— 30 v magnetic field— 16 gauss (5 amp) 

arc current— 1.5 amp flowrate— 1.6 gm/hr 


The above energy distribution is not an unreasonable 
distribution to assume for the electrons in the plasma. 
The primary electrons are those electrons from the cath- 
ode which are accelerated through the cathode sheath 
and undergo no collisions. The magnetic field of the 
engine will tend to randomly orient the directions of the 
velocity of the primary electrons. Maxwell electrons are 


2 Actually, one should analyze only the electron current to the probe. 
That is, use 

U = i + | U | 

where i e is the electron current, f is the probe current, and u is 
the saturated ion current. Since | U | < < * for most of the voltage 
ranges considered, the probe current has always been used instead 
of the electron current for the analysis. 


produced from electron-electron collisions between the 
primary electrons which have undergone inelastic col- 
lisions with atoms and the electrons ejected from the 
atoms by the ionization process. 



Fig. II. Corrected and uncorrected semilogarithmic LPC 
(operating conditions — same as for Fig. 10) 


Electron-electron and electron-atom elastic collisions 
will tend to produce a spread in the energy of the pri- 
mary distribution. One would expect the primary energy 
distribution to be better represented as gaussian or half- 
gaussian than as the monoenergetic spike used above. 
Thus, the current-voltage probe curves would not be 
exactly straight lines at large negative voltages, but a 
more complicated function, approximating a straight line 
as the gaussian becomes narrower. 

All of the data taken were analyzed by using the 
monoenergetic-plus-Maxwell energy distribution approxi- 
mation, even though in some cases the current-voltage 
curves were far from linear at negative probe voltages. 
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It must be emphasized that the assumed distribution is 
only an approximation of the actual condition and is used 
at this time because of its simplicity in application. A 
better method would be to compute the energy distribu- 


tion directly from the LPC as was done by Druyvesteyn 
(Ref. 6) and Medicus, but this method involves taking 
second derivatives of the characteristics, which is a 
highly inaccurate and tedious procedure. 


V. RESULTS 


Figure 12 is a drawing of the positions at which meas- 
urements were made inside the engine. The axial coordi- 
nates are labeled 1 through 5, while the radial coordinates 
are labeled A through E. 
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Fig. 12. Positions at which measurements were taken 

Figures 13 through 20 present the results of applying 
the foregoing analysis to the LPCs. These graphs show 
the axial and radial plasma density and energy variations 
under various operating conditions. The graphs appear 
in pairs, the first being the no-beam condition, and the 
second being the corresponding data when a beam was 
being extracted. Therefore, a comparison of static and 
operating conditions can be made. All of the plasma 
properties obtained from the probe curves have been 
plotted except the plasma potential, which was very 


difficult to measure. The plasma potential usually varied 
between 2 to 7 v above anode potential. The accuracy of 
these curves is estimated at only about 30%, the large 
error being due mostly to the straight-line approximation 
to the primary electrons. Even so, the trends can be 
determined. 

It is seen that the density of primary electrons drops 
off almost linearly with radial distance, while the Max- 
well density usually drops much more slowly. There is 
also a definite drop in electron temperature in the radial 
direction. The drop in the electron temperature may be 
explained if there is significant energy exchange between 
the electrons of the primary and Maxwell groups. Since 
the primary and Maxwell densities are larger near the 
center of the engine, the greatest energy exchange be- 
tween the groups would occur there; and since the energy 
of the primary electrons is higher than the Maxwell, the 
Maxwell mean energy would be raised. 

The number density of ions was not measured, but if 
the plasma is assumed electrically neutral, the number of 
ions should equal the total number of electrons. This 
number also falls off with radial distance, most often 
following the density of Maxwell electrons, since the 
ratio of Maxwell to primary electrons is usually low. 
In fact, in a few cases (Fig. 15) the density of Maxwell 
electrons was found to increase going out from the cath- 
ode, approach a maximum about halfway between anode 
and cathode, and then fall off quite rapidly. The total 
number of electrons (and therefore the ions) were concen- 
trated in .an annular ring about the cathode. An annular 
ring of increased intensity could sometimes be observed 
in the blue glow from the end of the engine. 
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N | = 0 15 x 1 0 M / cm 3 - NUMBER OF PRIMARIES AT A 
A/ z * 1.7 x IO M /cm 3 = NUMBER OF MAXWELLS AT A 
/V 3 = 1.85 x IO M /cm 3 - TOTAL NUMBER AT A 
N m -- NUMBER OF MAXWELL ELECTRONS 
N p = NUMBER OF PRIMARY ELECTRONS 
N t = TOTAL NUMBER OF ELECTRONS 




RADIAL POSITION 


RADIAL, SLOT POSITION 


RADIAL POSITION 


Fig. 13. Radial plot of electron density and energy variations (no accelerating potential applied) 

arc voltage — 50 v magnetic field - 20 gauss (6 amp) 

arc current— 4 amp flowrate— 1.5 gm/hr 
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/V, = 2.7 x I0 l0 /cm 3 = NUMBER OF PRIMARIES AT A 
N z = 3.9 x IO ,0 /cm 3 = NUMBER OF MAXWELLS AT A 
N z = 6.6 x lo'/'cm 3 = TOTAL NUMBER AT A 
N m = NUMBER OF MAXWELL ELECTRONS 
N p = NUMBER OF PRIMARY ELECTRONS 
N f = TOTAL NUMBER OF ELECTRONS 



RADIAL POSITION 


RADIAL, SLOT POSITION 


RADIAL POSITION 


Fig. 14. Radial plot of electron density and energy variations (accelerating potential applied) 

arc voltage— 50 v acceleration voltage— 5 kv 

arc current— 4 amp beam current (engine)— 160 ma 

magnetic field— 20 gauss (6 amp) bias voltage-1 kv 

flowrate— 1.5 gm/hr interception current— 0.8 ma 
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Fig. 15. Radial plot of electron density and energy variations (no accelerating potential applied) 

arc voltage-27 v magnetic field-17 gauss (5.2 amp) 

arc current-0.9 amp flowrate-1.4 gm/hr 
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yV, = LI 5 x IO ,0 /cm 3 = NUMBER OF PRIMARIES AT A 
N z = 0.409 x IO lo /cm 3 = NUMBER OF MAXWELLS AT A 
/V, = 1.56 x »O lo /cm 3 s TOTAL NUMBER AT A 
N m = NUMBER OF MAXWELL ELECTRONS 
N p = NUMBER OF PRIMARY ELECTRONS 




RADIAL POSITION RADI AL, SLOT POSITION RADIAL POSITION 

Fig. 16. Radial plot of electron density and energy variations (accelerating potential applied) 

arc voltage— 28 v acceleration voltage 4.9 kv 

arc current— 0.85 amp beam current (engine)— 32 ma 

magnetic fieid-1 7 gauss (5.2 amp) bias voltage— 1 kv 

flowrate— T. 4 gm/hr interception current— 0.4 ma 
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/V, = 0.508 x IO lo /cm 3 ,= NUMBER OF PRIMARIES AT A 
N 2 = 2.52 x IO ,0 /cn/ = NUMBER OF MAXWELLS AT A 
*3 = 3.03 x lO ,0 /cm 3 = TOTAL NUMBER AT A 
N m - NUMBER OF MAXWELL ELECTRONS 
N fi = NUMBER OF PRIMARY ELECTRONS 
N t - TOTAL NUMBER OF ELECTRONS 




RADIAL, GRIDS 


RADIAL POSITION 


Fig. 18. Radial plot of electron energy and density variations near grids (accelerating potential applied) 

(operating conditions — same as for Fig. 14) 
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A/, = 1.09 x IO ,0 /cm 3 = NUMBER OF PRIMARIES AT 5 
/V 2 = 3.7 x I0'%m 3 = NUMBER OF MAXWELLS AT 5 
/V 3 = 4.8 x IO f0 /cm 3 = TOTAL NUMBER AT 5 
N m = NUMBER OF MAXWELL ELECTRONS 
N p = NUMBER OF PRIMARY ELECTRONS 
N t = TOTAL NUMBER OF ELECTRONS 



AXIAL POSITION AXIAL, POSITION C AXIAL POSITION 

Fig. 19. Axial plot of electron energy and density variations (no accelerating potential applied) 

arc voltage— 30 v magnetic field-16 gauss (5 amp) 

arc current— 1 .5 amp flowrate— 1.4 gm/hr 
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2 0 




yv, = 1.17 x IO ,0 /c™ 3 = NUMBER OF PRIMARIES AT A 
N 2 = 1.5 x IO ,0 /cm 3 = NUMBER OF MAXWELLS AT A 

yv 3 = 2.67 x IO ,0 /crn 3 = TOTAL NUMBER AT A 
JV m = NUMBER OF MAXWELL ELECTRONS 
N p - NUMBER OF PRIMARY ELECTRONS 
Nt- TOTAL NUMBER OF ELECTRONS 




Fig. 20. Axial plot of electron energy and density variations (accelerating potential applied) 


arc voltage— 30 v 
arc current— 0.9 amp 
magnetic field— 16 gauss (5 amp) 
flowrate— 1.4 gm/hr 


acceleration voltage— 4 kv 
beam current (engine)— 40 ma 
bias voltage— 1 kv 
interception current— 0.4 ma 
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The effect of applying accelerating voltage is seen to 
reduce the density of Maxwell electrons (at the particular 
transverse position at which the radial traces were taken), 
which drastically increases the ratio of primary to Max- 
well electrons. It is not clear why the accelerating poten- 
tial should cause such an effect. One would expect that 
the plasma near the grids would shield the interior 
plasma quite effectively, but evidently this is not the case. 
The extraction of ions through the grids must set up axial 
electric fields which alter the spatial distribution of the 
lower-energy Maxwell electrons. The axial spatial distri- 
butions taken from axial sweeps of the Langmuir probe 
(Fig. 19, 20) do in fact show a gross rearrangement of the 
Maxwell electrons when high voltage is applied. 


To see the effect of the primary electrons upon plasma 
formation, consider 


N + (r, E ) 


dN + (r, E) 
dt 


the rate of production of positive ions at any point r by 
electrons of energy E. This is given by 

$+ ( r >E) = N 0 (r)f(E, r)°i(E)v (5) 

where 

N a (r) = the number density of neutral atoms at the 
point r 

f(E, r) — the number of electrons with energy E at 
the point r, i.e., the spatial energy distribu- 
tion of the electrons 

<^(E) = the cross section for ionization of a neutral 
atom 

v = the velocity 7 of the electron 

Assume that the energy distribution of electrons is made 
up of two parts, a primary and a Maxwell part, so that 
f can be written: 

m r) = UE) N„ (r) + f m (E)N m (r) (6) 

where the subscript p stands for primary and m for 
Maxwell; N p and N m are the spatial number density 
distributions; and f p and f m are the energy distributions. 
For the primary energy distribution, assume a mono- 
energetic distribution at an energy E p . Then, 

f p (E) = S(Ep) = Dirac $ function (6') 

The rate of ion production at point r is determined by 
integrating (5) over all energies, using the energy distri- 
bution (6), (6') and the usual Maxwell energy distribution: 


N + (r) 


N 0 (r) [f t (E)N p (r) + f m (E)NJr)] 
°,(E)\fEdE 


This yields on integration 

N + (r) = N,N 0 (r)N P „(r)VE; [(? p (r) + Q„(r)] 

( 7 ) 

where 

77i — electron mass 

o m — maximum value of the electron-mercury 
ionization cross section 

N x - Np/N pn 

= maximum value of primary density 
N pn (r) = normalized primary density 

E p ~ primary energy 
Q p (r) = ° n (E„(r)/V % ) 

— normalized ionization cross section of mer- 
cury; the approximation of ° n used was 
developed by Kerrisk (Ref. 2 and 7) and is 
shown in Fig. 21. 



Fig. 21. Approximation of the mercury ionization 
cross section 
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ionization coefficient of Maxwell elec- 
trons; this function has been evaluated 
by Kerrisk and is reproduced in Fig. 
22 

mean Maxwell energy 

first ionization potential of mercury 



Fig. 22. Ionization coefficient for 
Maxwell electrons 


Maxwell electrons, the function Q m (r) is at most Mo as 
large as the function Q p (r). Furthermore, Q m is a sensi- 
tive function of the Maxwell mean energy and, since the 
mean energy always falls off in the radial direction, the 
ratio of Q m to Q v will be smaller at the anode than it is 
at the cathode. In other words, under most conditions, 
the term Q m in Eq. (7) can be neglected in comparison to 
Q p . It follows that the primary electrons are more effec- 
tive in producing ions. 



RADIAL POSITION 

Fig. 23. Ionization rate curve (no accelerating 
potential applied) 

(operating conditions — same as for Fig. 13) 

The ion current density J + necessary to form a stable 
sheath is given by Bohm (Ref. 8) to be 

’• = eN - VW < s > 



Note that the functions Q p and P m will only have a 
radial variation if E p and E m vary radially. The functions 
Q p and Q m of Eq. (7) have been calculated from the 
spatial density distribution curves of Fig. 13 and 14 and 
are shown plotted in Fig. 23 and 24. One sees immedi- 
ately that even though there are fewer primary than 


If it is assumed that this equation is valid in the extrac- 
tion region of the engine near the grids, then from the 
measured values of the total ion density ( N + ) and mean 
Maxwell energy (T e ), the ion current can be calculated. 
The calculation has been performed on the probe data 
obtained when the engine was producing a 160-ma beam 
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(Fig. 17, 18) with agreement to the corrrect order of 
magnitude. At the grids the current density near the 
anode was about Vr the current density near the center 
of the engine. 

All of the data presented in this section are intended 
to be representative of a large amount of data under 
various operating conditions obtained from the probe 
curves. All the probe data for the following variations 
in engine parameters yield essentially similar results for 
the shape of the spatial density distributions. 


Arc voltage 
Arc current 
Magnetic field 
Cathode temperature 
Propellant flowrate 
Net acceleration voltage 
Bias voltage 
Beam current 


20 to 60 v 
0.85 to 4.7 amp 
13 to 23 gauss 
950 to 1250 °C 
1.65 to 2.2 gm/hr 
3 to 8 kv 
1 to 1.5 kv 
40 to 300 ma 



RADIAL POSITION 

Fig. 24. Ionization rate curve (accelerating 
potential applied) 

(operating conditions — same as for Fig. 14) 
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VI. CONCLUSIONS 


Langmuir probe curves analyzed for the electron 
energy distribution for an operating (30 to 160-ma beam) 
and static (no beam, but plasma) Kaufman-type ion 
engine show that the energy distribution is not entirely 
Maxwellian. The energy distribution was found to con- 
sist of a high-energy primary distribution and a lower- 


energy Maxwell distribution. The primary electrons 
appear to be the main source of ion production. Further, 
the ion density near the grids as well as in the interior 
of the engine is not uniform. Any theoretical model of 
this source must include the primary distribution of 
electrons. 
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